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We investigate the scale-dependence, or the runnings, of linear and second order den- 
sity perturbations generated in various curvaton scenarios. We argue that the second order 
perturbations, i.e. non-Gaussianity, can strongly depend on the scale, even when the linear 
perturbations are nearly scale-invariant. We present analytic formulae for the runnings from 
curvatons with general energy potentials, and clarify the conditions under which /nl becomes 
strongly scale-dependent. From the point of view of the /nl running, curvaton potentials can 
.£h be classified into roughly two categories by whether the potential flattens or steepens com- 

pared to a quadratic one. As such examples, we study pseudo-Nambu-Goldstone curvatons, 
and self-interacting curvatons, respectively. The dynamics of non-quadratic curvatons and 
the behaviors of the resulting density perturbations are clarified by analytical methods. Then 
we also study models where multiple source can be responsible for density perturbations such 
as the multi-curvaton, and mixed curvaton and inflaton models where the running of /nl 
can also be large due to their multi-source nature. We make quantitative analysis for each 
curvaton scenario and discuss in what cases the scale-dependence, in particular, of /nl can 
be large enough to be probed with future CMB experiments. 
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1 Introduction 

The origin of density fluctuations in the Universe is one of the important issues in cosmology 
and also gives invaluable information for high energy physics since they are assumed to be 
generated in the very early Universe. Although quantum fluctuations of the inflaton, which 
drives inflation in the very early Universe, has been considered to be its origin over the 
years, other mechanisms have also been discussed. Among possible candidates, the curvaton 
scenario [1-3] has been attracting much attention due to several reasons, one of which is 
the recent observational results on primordial non-Gaussianity. How much the primordial 
fluctuations deviate from a Gaussian distribution can be characterized by the non-linearity 
parameter / NL , whose constraint from current observations is —10 < /nl < 74 (95 % C.L.) [4] 
for the so-called local type non-Gaussianity* 1 . If future data such as those from the on-going 
Planck satellite confirms large local type non-Gaussianity of primordial fluctuations at the 
level of /nl > ^(10), it readily excludes conventional single-field inflation models predicting 
/nl *C Oil) as the origin of the density perturbations. On the other hand, the curvaton can 
generate large non-Gaussianity at this level, which would motivate us to seriously consider this 
scenario. Furthermore, the curvaton may be well fitted naturally into some particle physics 



^ tl In this paper we mainly discuss local type non-Gaussianities. For constraints on /nl of other types, 
see [4]. 



or string theory models (see, for e.g., [5-17]), thus also in this respect, the curvaton has been 
the target of intense study. 

In many works on the curvaton model, it is assumed that the potential for the curvaton has 
a simple quadratic form and the curvaton is totally responsible for cosmic density fluctuations. 
However, such assumptions are removed in some microscopic constructions of the curvaton. 
Regarding the former assumption, there exist models realizing potentials that deviate from a 
purely quadratic one. The consequences of such non-quadratic curvaton potentials have been 
investigated for the self-interacting curvaton [18-24] and the pseudo-Nambu-Goldstone (NG) 
boson one [25-28] . Interestingly, it has been shown that the predictions for non-Gaussianity of 
the primordial perturbations significantly depend on the form of the potential. In particular, 
when the potential deviates from the quadratic form, the non-linearity parameter /nl (and 
also other parameters characterizing non-Gaussianity such as r NL and #nl for the trispectrum) 
can be considerably scale-dependent [29-33], even when the linear order perturbations are 
nearly scale- invariant. On the other hand, the assumption of the curvaton being totally 
responsible for the density perturbations is removed for multi-source scenarios, where the scale 
dependence of non-Gaussianity can arise even with quadratic curvaton potentials [29,30]. Such 
a situation can be realized in the mixed curvaton and inflaton model [34-38] where fluctuations 
from the curvaton and inflaton can be both responsible for cosmic density perturbation today. 
One could also consider a model where there exist multiple curvaton fields and they can be 
responsible for density fluctuations. Some authors have studied a model of two curvatons 
[39,40]. 

It should also be noted that in the event of a detection of large non-Gaussianity, upcoming 
CMB experiments can set severe constraints on the scale-dependence of /nl as well, especially 
when combined with large-scale structure surveys, see e.g. [41-44]. Such information beyond 
/nl would be a powerful probe of the physics of the early universe. In the light of these consid- 
erations, the issue of the scale-dependence or "running" of non-Gaussianity is an interesting 
and important subject. 

In this paper we carry out analytical studies of density perturbations from curvatons, 
by applying the method developed in [28]. We derive generic formulae for the runnings of 
the linear and second order density perturbations, which enable us to go beyond individual 
case studies and give a systematic treatment of the scale-dependence of /nl in the curva- 
ton mechanism. Our study not only serves as an analytical counterpart to previous works 
that basically relied on numerical computations, but also clarifies the underlying reason why 
strongly scale-dependent non-Gaussianities can be generated from non-quadratic curvatons. 
Furthermore, we present conditions for curvatons to produce large running of /nl, which 
turn out to take simple forms (2.34) and (2.35) for a single curvaton. In passing, we also 
investigate the running of the spectral index of the power spectrum a, which can also give 
interesting signatures. We will clarify the direct relation between a and the running of /nl, 
which we parametrize as n / NL . The general study of n^ NL and a would give new insight into 
the curvaton mechanism. 



Curvaton potentials can roughly be categorized by whether they flatten or steepen com- 
pared to a quadratic one. It will be shown that this classification is important for discussing 
rif NL from curvatons. As typical examples of the two cases, we will look into pseudo-Nambu- 
Goldstone curvatons and self-interacting curvatons, respectively. For the former case, the 
relation between rif NL and a plays an important role, restricting the running of /nl from 
current observational bounds on a. On the other hand, for the latter case, strongly scale- 
dependent / NL can be produced even for a suppressed a due to the steepness of the potential. 
Regarding n/ NL generated in multi-source models, we make a quantitative study for some 
models and show in what cases the scale dependence can be large enough to be probed in 
future cosmological observations. 

The structure of this paper is as follows: In Section 2 we derive analytic expressions for 
the scale-dependence of the linear and second order density perturbations from a curvaton. 
Then we move on and apply the generic discussions to pseudo-Nambu-Goldstone curvatons 
in Section 3, and self-interacting curvatons in Section 4. In Section 5, we further investigate 
the issue for the mixed curvaton and inflaton scenario. Then in Section 6, we discuss multi- 
curvaton model. We present our conclusions in Section 7. 

We give a brief discussion on density perturbations from curvatons with non-sinusoidal 
oscillations in Appendix A. 

2 Scale-Dependence of /nl m the Curvaton Mechanism 

A light curvaton acquires nearly scale-invariant field fluctuations during inflation, that are 
converted into the cosmological density perturbations as the curvaton oscillates and decays in 
the post- inflationary era. In this section we derive generic expressions for the running of non- 
Gaussianity from a curvaton, and discuss the conditions under which the curvaton sources a 
largely scale-dependent /nl- Our discussions in this section are based on the work [28] which 
developed analytic methods for computing density perturbations in the curvaton mechanism. 
We extend their results and compute the runnings of the linear and second order perturba- 
tions. In this Section and also in the following Sections 3 and 4, we assume that density 
perturbations sourced from the inflaton are neglected. However, we remove this assumption 
in Section 5 where we study the mixed curvaton and inflaton scenario. 

2.1 Density Perturbations 

Let us start by laying out some results of [28] for density perturbations produced by a cur- 
vaton a possessing an effective potential V(a). The potential is assumed to have no explicit 
dependence on time, and also that it is well approximated by a quadratic one around its 
minimum so that the curvaton oscillations are sinusoidal* 2 . The curvaton energy density is 



# 2 The formulae can be generalized to cases with non-sinusoidal oscillations as well, see Appendix B of [28] 



considered to redshift similarly to nonrelativistic matter after the onset of the oscillations 
until when the curvaton decays into radiation, whereas we suppose the infiaton to behave as 
matter from the end of inflation until reheating (=inflaton decay), after which the infiaton 
turns into radiation. The energy density of the curvaton before the beginning of its oscilla- 
tion is assumed to be negligibly tiny compared to the total energy of the Universe, having 
little effect on the expansion history. Furthermore, the Hubble parameter during inflation is 
considered to be nearly constant. For detailed discussions on the derivations of the following 
results, we refer the reader to [28]. 

Using the 5Af- formalism [45-48], the density perturbations are obtained by computing 
quantities such as the time when the curvaton oscillation starts and the curvaton energy 
density as functions of the curvaton field value <r* (hereafter the subscript * denotes values 
when the CMB scale fc* exits the horizon). The curvaton dynamics prior to the oscillation 
can be tracked by the attractor solution 

{3 (during inflation) 
9/2 (matter domination) (2.1) 

5 (radiation domination) , 

which is a good approximation while \V" /cH 2 \ <C 1. Here, a prime denotes a derivative with 
respect to a, and an overdot a time derivative. Setting the minimum of the potential about 
which the curvaton oscillates to a = 0, we can define the onset of the oscillation as when the 
time scale of the curvaton rolling becomes comparable to the Hubble time, i.e. \a/Ha\ = 1. 
This, combined with (2.1), gives the Hubble parameter at the time, 

HL = ^4 (2.2) 

where the subscript "osc" denotes values at the onset of the curvaton oscillation, and c is a 
constant whose value is set by whether reheating (= infiaton decay, at t re h) is earlier/later 
than the onset of the curvaton oscillation: 
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We define the power spectrum V^ of the density perturbations ( as 

271-2 

(Ch(k>) = {2'K) 3 6®{k + k , )Pc{k) with P c (k) = — V c (k), (2.4) 

where k = \k\. Supposing the curvaton field fluctuations to satisfy V$ a (k) = (H \k=aH/^) 2 
at the time when the scale k exits the horizon, then the linear order density perturbations at 
the CMB scale can be expressed in terms of the curvaton potential as 



with 

dH _ r M w vwi/ ngonj 3X(cr osc ) | V> osc ) ,- R . 
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Here, f is the energy density ratio between the curvaton and radiation (which originates from 

the inflaton) upon curvaton decay* 3 
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(2.8) 

dec 



while the function X denotes effects due to the non-uniform onset of the curvaton oscillations 
(which are absent for a purely quadratic curvaton potential), defined as follows: 

*<»-) = ^3) { y , (0 - Ij • (2-9) 

From the bispectrum 

(CfciCfcaCfc.) = (27r)V 3) (fc 1 + fc 2 + fc 3 ) J B c (A; 1 ,A ;2 ,A;3), (2.10) 

one can generally define the non-linearity parameter /nl as 

S c (fci, k 2 , h) = ^ NL (fc l5 fc 2 , h) [P c (A;i)P c (A; 2 ) + P c (k 2 )P ( (k 3 ) + P^k^P^k,)] . (2.11) 

However, upon parameterizing the overall amplitude of the local-type bispectra produced by 
curvatons,* 4 throughout this paper we discuss /nl on the equilateral configuration 
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f„ L (k) = -B ( (k,k,k) l—V ( (k)) . (2.13) 



We suppose that the bispectrum of the curvaton fluctuations (£cr fcl 5cr fc2 5cr fe g) with |fc 1 |=|fc 2 |= 
|fe 3 | = A; vanishes when k exits the horizon, and further neglect direct trispectra of 8a (i.e. 
(5a qi 5(Tq 2 5a qa 5a q4 ) = (5a qi 5a q2 )(5a q3 5a q4 ) + (2 perm.)). Then it can be shown that /nl 



# 3 The following quantity is also used in the literature to express the formulae for density perturbations from 
the curvaton: 

3yO CT 
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4p r + 3p a 



(2.7) 



which is also evaluated at the curvaton decay. Notice that f defined here is a bit different from fdec- 

# 4 It should also be noted that, especially when /nl is strongly scale-dependent, the bispectrum (2.10) 
sourced from a curvaton has shapes similar to, but not exactly of the "local form" [49] which is often given as 

Si oca i(fc 1 ,fc 2 ,fc 3 )cxp-3+p-3+p-3. (2.12) 



takes the form [28] 
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The scale-dependence of /nl shows up through cr*, whereas <t osc and f are independent of 
the wave number. Hence for later convenience, let us divide the expression into terms that 
explicitly depend on a*, and the rest, 



/nl = fi{<r*, ow, r) + / 2 (o"osc, r), 



where we have introduced 
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Upon obtaining the second line, we have used (2.6)* 5 . 
# 5 For a curvaton with a quadratic potential V oc cr 2 , one obtains 
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Furthermore, the terms /i and /2 depend only on f, 

5(4 + 3f) 
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/a = 
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and thus 
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Large non-Gaussianity is generated for f<l, under which 

/nl 



f - 1 . ~ 5 
/l -2 /2 "37- 



(2.15) 



(2.16) 



(2.17) 

(2.18) 
(2.19) 

(2.20) 



The energy density ratio f (2.8) is obtained as 

f = Max 
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where the first and second terms in the Max. parentheses correspond to the curvaton being 
subdominant and dominant at its decay, respectively, while the Min. parentheses are due 
to whether the onset of oscillation is after or before reheating. 1^ and IV are constants 
denoting respectively the decay rates of the inflaton and the curvaton. Throughout we adopt 
the sudden decay approximation where the scalar fields suddenly decay into radiation when 
H = T. 

The curvaton field value at the onset of the oscillations <r osc is obtained by integrating (2.1), 

d ° K 1 (2.22) 



V 3HI 2c(c-3)tf 2 sc ' 

solving which gives <t osc as a function of a** 6 . Here, A/"* is the number of e-folds during 
inflation between the horizon exit of the CMB scale and the end of inflation, and H m { is 
the inflationary Hubble scale (we are assuming a nearly constant Hubble parameter during 
inflation, thus H{ n { ~ H*). Let us also show the derivative of <7 0SC , 

^ = (l-X(Or 1 ^, (2-23) 

which has entered (2.6) through dN /da* oc d<T osc /da*. 

Thus by combining the above expressions, one can compute the resulting density pertur- 
bations from a curvaton with a generic potential V(a), given the curvaton field value at the 
CMB scale horizon exit a*, the decay rates of the inflaton T^ and curvaton r CT , the inflationary 
scale H iQ {, and the duration of inflation A/"*. 

2.2 Scale-Dependence 

Now we extend the above expressions to discuss the scale-dependence of the linear and second 
order density perturbations, which are the main topic of this paper. 

Since we are assuming \H/H 2 \ ^C 1 during inflation, the comoving wave number k at 
around the CMB scale satisfies 

d\nk~H*dt. (2.24) 



* 6 When (2.22) admits as solutions for tr osc both positive and negative values, one should take the sign of 



cr osc to match with that of a* . 



Then by using the slow-roll approximation for the curvaton 

3if„<7, ~ -V'fa), (2.25) 

one obtains the spectral index of the linear order perturbations at the CMB scale 



-— <L»«° 2 fH^ 



as well as its running 



dink H* Ht 3H% HI 9 Hi l ' ' 

Focusing on the contributions to the scale-dependence that are sourced purely by the tilt of 
the curvaton potential, we introduce the following parameters 

& = ~g ^ J H 4 - • (2-29) 

When the Hubble parameter during inflation is exactly a constant, these give the spectral 
index and its running at the leading order, i.e. n s ~ n s and a ~ a. 

In order to parametrize the scale-dependence of the non-Gaussianity, we define the spectral 

index of /nl as follows: 

_ d\n\f NL \ 

n *» = -dwT- (2 ' 30) 

We remark that a totally scale- invariant /nl corresponds to n/ NL = (instead of 1). Then, 
in a similar fashion as above, one arrives at 

1 5(4 + 3r) f V'(a osc ) 3X(a 0SC ) \-\ V'(* m ) V'{tx m ) 

"/nl - 1 T^ \ 777 7 f C 1 _ -M^oscJJ T777 r — ^— • (2.31) 

/nl 18r I V(a osc ) a osc J V"(ct osc ) Hi 

Here we note that time derivatives of H* do not show up at the leading order, since /nl (2.14) 
does not explicitly depend on the Hubble parameter during inflation. One clearly sees that a 
no n- vanishing n/ NL (at the leading order) requires non-zero V'"(a*), and hence /nl produced 
by a curvaton with a quadratic potential is scale-invariant. This opens up the possibility 
that a slight deviation of the curvaton potential from a quadratic one can be verified through 
observing the running of /nl- 

By using (2.6), one can simplify (2.31) as 



which shows that when the power spectrum is fixed to a certain value, e.g. from the COBE 
(WMAP) normalization, then the value of the product n/ NL /NL is determined only by infor- 
mation at the CMB scale horizon exit. 

Moreover, in terms of the parameters (2.16), (2.28), and (2.29), the expression (2.31) can 
be recast into the form of 

nfia. - ^T A- (2-33) 

n s — i Jnl 

This equation makes clear the typical amplitude of the running of /nl expected in the curvaton 

model. In simple cases, the running parameter a (though not necessarily equivalent to the 

actual running of the spectral index) has a smaller size than the spectral index parameter n s — 

1. Furthermore, since f\ explicitly depends on a* while the rest of the terms of /nl do not* 7 , 

one may not expect /i to be much larger than the sum /nl = /i + /2- Thus naively one would 

expect the running of /nl from curvatons to be highly suppressed. 

In other words, nf NL as large as, say |n/ NL | > 1, is realized only when at least one of the 

following two conditions are satisfied: 

|a|>|n.-l|, (2.34) 

IM>I/nl|- (2.35) 

The former condition (2.34) may be realized by curvaton potentials possessing inflection 
points, or more generally, by potentials that flatten compared to a quadratic as one goes away 
from the minimum. As such an example, in Section 3 we study pseudo-Nambu-Goldstone 
curvatons with cosine-type potentials. However, there we will see that the resulting rif NL is 
actually directly bounded by observational constraints on the running a of the spectral index 
of the linear order perturbations. The latter condition (2.35) indicates that f\ and f'2 cancel 
each other, suppressing /nl compared to /1. Such suppression of /nl is known to exist for 
self-interacting curvatons possessing polynomial terms that are higher order than quadratic, 
as was numerically shown in [18-24]. In Section 4 we look into self-interacting curvatons and 
see that large n/ NL can be obtained even under a suppressed running a, and further discuss 
that such behavior stems from the curvaton potential steepening more rapidly than quadratic 
ones. 

The value of /nl itself needs to be large for its running to be detectable, and thus experi- 
ments are sensitive to the product ?t./ nl /nl, instead of rif NL alone. Upon discussing example 
models in the following sections, we refer to the results of [42] where detectability of a scale- 
dependent non-Gaussianity through CMB experiments are analyzed* 8 . The running of the 

# 7 Strictly speaking, <7 OS c also depends on cr* through (2.22), but the point here is that f\ and fa depend 
quite differently on a. t . 

# 8 We note that the definition of our running parameter n/ NL (2.30) is not exactly the same as the n^c 
parameter adopted in [42] , especially in the sense that ting is set to a constant whereas Uf Nh itself can run. 
Nevertheless we adopt their results upon discussing detectability. 

9 



local-type /nl can be probed if it is large enough to satisfy 

50 
|n /NL |>Ax— , (2.36) 

Jnl 

where the right hand side takes A ~ 0.68, 0.10, and 0.05 for WMAP [4], Planck [50], and 
CMBPol [51], respectively, assuming as fiducial values /nl = 50, rif NL = 0, and a full-sky 
coverage. 

Before ending this section, we should mention about errors in the analytic formulae. 
Firstly, the approximation (2.24) contains error of order H*/H 2 , and (2.25) of order V"(cr*)/H 2 . 
By taking into account such errors, one can check that they do not modify the above results 
onfjj, a, and rif NL at the leading order. 

However further approximations and simplifications have been carried out upon obtain- 
ing the analytic expressions of Subsection 2.1, e.g., the approximation (2.1) on the curvaton 
dynamics, which is correct up to order V" /H 2 . This can source errors in the results of order 
V" /H 2 , and also of derivatives of V" /H 2 in terms of cr* and/or time t. It should be noted 
that such error with various orders of derivatives can accumulate and lead to breakdown of 
the analytic expressions especially for higher-order correlation functions and running. Nev- 
ertheless, for the explicit examples we study in the following sections, we will see that the 
analytic results match well with results from numerical computations. 

3 Pseudo-Nambu-Goldstone Curvatons 

Since the parameters n s and a in the expression for rif NL (2.33) are not necessarily the actual 
spectral index and its running under a non-vanishing H during inflation, one may expect that 
the condition (2.34) for a large rif NL is easily satisfied without contradicting with observa- 
tional constraints on the flatness of the power spectrum. One may imagine cases where the 
amplitudes of a and n s — 1 are much smaller than unity, though possessing a hierarchy among 
them as \a\ 3> \n s — 1|. 

In this section we study the case where the curvaton is realized as a pseudo-Nambu- 
Goldstone (NG) boson of a broken U(l) symmetry, possessing a cosine-type potential [25-28]. 
Given that the curvaton at the CMB scale horizon exit is located close to the inflection point 
of the potential, h s — 1 vanishes while a remains finite. However, we will see that the resulting 
rif NL is actually set by the absolute value of a, thus bounded by observational constraints on 
the running of the spectral index (unless a is cancelled out by the H terms in (2.27)). 



We consider the potential of the form 



V{a) = A 



I - cos I - 



10 



(3.1) 



where / and A are mass scales. It has inflection points at a j f = (1/2 + n)ir with n G Z, and 
we focus on the region < a osc / f, a*/ f < n without loss of generality. Since 

h '-^lwp cos {j)- (3 ' 2) 

. 2 A 8 . 2 (<j,\ 

one sees that a/(n s — 1) blows up at er*// = ir/2. However, n/ NL is actually set by the value 
itself of a as we will soon see. 

In Figures 2 - 7, we display the density perturbations plotted as a function of er*// under 
the parameter set A = 10 15 GeV and / = 10 17 GeV, along with the inflationary parameters 
Hi n f = 3.5 x 10 13 GeV, T^ = 10 11 GeV (i.e. the energy density at reheating (= inflaton 
decay) is p reh ~ 6.5 x 10 14 GeV), and A/"* = 50. # 9 The curvaton decay rate is set to T a = 
Y^ — U — — i^ts") supposing that the coupling of the NG curvaton with its decay product 
is suppressed by the symmetry breaking scale /. This set of parameters are chosen such that 
the COBE (WMAP) normalization value [4] V^ ~ 2.4 x 1CT 9 as well as n s ps 0.96 are realized 
at around o^jnj w O.8.* 10 Moreover, the curvaton starts its oscillation before reheating, thus 
c = 9/2. 

cr osc is computed by solving (2.22), which now takes the form 



In 



tan (cr osc /2/) 
tan(a*/2/) 



AC A 4 1 a osc /f 



3Hl f P 2(c-3)sin(a osc //)' 



By obtaining a fitting function for the solution <j osc (<j*) (which is shown in Figure 1), we have 
analytically calculated the density perturbations in terms of cr*. The analytically estimated 
results are shown as blue solid lines in Figures 2-7 for the region 0.01 < cF*/^f < 0.99. 

We have also numerically computed the density perturbations, by solving the curvaton's 
equation of motion and computing the differences in the number of e-foldings obtained from 
different initial values a*. Upon carrying out the numerical computations, we have set the 
inflaton energy to a constant during inflation which lasts for 50 e-foldings after the CMB scale 
exits the horizon, then transferred the inflaton energy to non-relativistic matter redshifting 
as p oc or 3 , and finally to radiation as p oc a~ 4 . Furthermore, we adopted the sudden decay 
approximation for the inflaton/curvaton, as well as a sudden end of inflation. The numerically 



# 9 The e-folding number AC is basically determined by knowing the scale of inflation and the subsequent ex- 
pansion history. However in order to clarify the dependence of the density perturbations on each parameter, 
we fix A/"* to 50 in Sections 3 and 4. The direct consequence of a larger (smaller) A/"* is to decrease (in- 
crease) Cose j hence docs not affect n s nor a, while other cosmological obscrvablcs can be affected, especially 
since f is also changed. Nevertheless, the overall behavior of the density perturbations are not influenced 
much by the detailed value of Af* ■ Thus the fixing of AC does not affect our discussions greatly. 
# 10 The inflationary scale needs to be rather high when one tries to realize the COBE (WMAP) normalization 
and \n s — 1| ~ 0.01 at c* values not so close to the hilltop [28]. For such high-scale inflation, depending on 
the inflationary mechanism, density perturbations from the inflaton can also become substantial. 
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computed results are shown as blue dots in the figures. Since we have fixed the inflaton energy 
density to a constant during inflation, the spectral index (2.26) and its running (2.27) are 
expected to match with the parameters n s (2.28) and a (2.29), respectively, at the leading 
order. Thus we have displayed n s and a, and the numerically computed n s and a in the same 
plots in Figures 2 and 3. One sees that the analytic estimations are in good agreement with 
the results of the numerical calculations in all Figures 2-7. 

In most of the displayed cr*// region the behavior of the density perturbations from a NG 
curvaton can be understood similarly as for a curvaton with a quadratic potential (cf. Foot- 
note #5), except for around the hilltop a* pa ixf where the linear perturbations and /nl are 
enhanced. (In the figures we have plotted up to a^/nf < 0.99, however as one goes even 
closer to the hilltop, V( and /nl further increase, cf. [28].) 

The amplitude of the running of the non-Gaussianity nj NL is more or less correlated with 
the at parameter, which is independent of the energy fraction f. This is a rather generic 
feature of density perturbations from a NG curvaton, which can be understood from (2.32), 

s i foAfy 1 a" . /v.\ / 6 aw v'», [& 

In the far right hand side, the product inside the parentheses is a combination whose amplitude 
is smaller than unity for quadratic curvatons.* 11 Here, recall that the density perturbations 
from a NG curvaton are more or less the same as those from a quadratic curvaton, except for 
the hilltop region. Therefore the () _1 term in the far right hand side of (3.5) is roughly of 
order unity (except for when /nl vanishes), and accordingly, |n/ NL | is roughly of order y/&. 
On the other hand, when |/nl| "C 1, the product inside the parentheses also becomes much 
smaller than unity and thus its inverse blows up (though not seen in Figure 6, this can 
happen depending of the parameter set). However, in such a case the product rif NL f NL to 
which experiments are sensitive is suppressed. Hence one can conclude that the running 
of /nl from a NG curvaton is basically set by the absolute value of a, which is constrained 
by current observational bounds on a running spectral index* 12 (unless the time- variation of 
the Hubble parameter cancels out a from (2.27)). 

As mentioned above, the produced perturbations behave quite differently as one ap- 
proaches the hilltop, i.e. cr*/ / — > it. However there the potential is well approximated 
by a quadratic, i.e. V ~ Vo — m 2 a 2 , hence at approaches zero and n/ NL is suppressed. 

We show contours in the n/ NL - /nl planes from a NG curvaton in Figures 8 and 9. Instead 
of choosing a certain parameter set for e.g. A and / as in the previous figures, here we have 



# n A quadratic curvaton V ex er 2 gives 

6 f dN 9f 2 -8f-16 , oc , 

-7/NL'n — -<r* = — — , (3.6) 

5 <9ct* (3r + 4:) A 

which monotonically increases with f, and approaches +(— )1 in the limit r 3> (<C)1. 
# 12 The current observational la limit for the running spectral index from WMAP7 is a — — 0.034±0.026 [4]. 
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Figure 1: cw as a function of a*. 



fixed n s to 0.96 and 1.04, respectively, and plotted the contours in the following way: For 
a given h s , by using (3.2) one can rewrite A 4 /H? n{ f 2 in terms of a*/ f. Hence /nl as well 
as rif NL are uniquely determined by a*//, cr 08c /f, f, n s , and c. Furthermore, a osc / f is given 
as a function of a*// after solving (3.4) (note that A A /H 2 n{ f 2 is now determined by n s and 
a *l f)i where we took A/"* = 50 and c = 9/2, i.e. t rc h > t osc (taking instead c = 5, i.e. 
^reh < ^oso niakes little difference for the resulting contours). Each colored solid line in the 
figures is plotted for a fixed value of cr^/nf, under a varying f . Different values for o^jixj give 
different values for a: blue, green, and red lines correspond to a ~ 0.05, 0.008, and 8 x 10~ 5 , 
respectively. We note that this procedure for plotting the contours does not fix the linear 
perturbation amplitude, i.e., the individual values for A, /, and H ini can be varied to change 
the perturbation amplitude for a fixed A 4 /H 2 nf f 2 . As reference values for detectability, we 
have also shown n/ NL /NL/50 = 0.10 contours as black dashed lines. Outside these contours 
corresponds to the region detectable by Planck, cf. (2.36). 

Taking smaller f (<C 1) corresponds to moving upwards along the colored contour lines. 
It is clearly seen that a large % nl /nl is realized for larger values of a, as was discussed below 
(3.5). Since the contours in the two figures are chosen to take the same values for a, they are 
more or less the same in the f <C 1 regime. 

On the other hand, as f is increased, the contour lines eventually turns left and rif NL 
decreases. There the blue contours realize |n/ NL | much larger than \a\, due to the suppressed 
|/nl| Cl In some cases /nl can even cross zero (as it happens for quadratic curvatons, 
cf. (3.6)), which is accompanied by nf NL blowing up and changing sign. This is seen in 
Figure 9 as the contour lines extending towards n/ NL — > — oo, and then coming back from the 
right. Although such behavior is absent in Figure 8, it can happen for a*// > tc/2 cases as 
well, depending on the parameter values. Since the blowing up of rif NL happens when /nl 
is tiny, its observational detection would be challenging, even if it happened. In the f — > oo 
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Figure 2: Curvaton contribution to the spec- 
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Figure 3: Curvaton contribution to the run- 
ning of the spectral index (2.29). 
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Figure 7: Running of non-Gaussianity. 
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limit, nf NL approaches a certain value, which can be seen as the end points of the contour 
lines. (The end points of the blue lines in both figures are outside the displayed region.) 

Focusing on a contour line with a fixed value for cr*/irf, then taking a larger value for 
\n s — 1| is equivalent to increasing a, hence the contour would shift towards larger |n/ N J. In 
summary, for NG curvatons, large |ti/ nl /nl| is produced together with a large a, therefore 
it is already strictly constrained by current observational bounds on running spectral index, 
unless a is cancelled out from the expression for a due to a varying Hubble parameter during 
inflation. 
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Figure 8: Varying f under h s = 0.96. cr*/nf 
is fixed to 0.55 (blue), 0.60 (green), and 0.90 
(red), corresponding toa~ 0.03, 0.008, and 
8 x 10~ 5 , respectively. The expected obser- 
vational sensitivity of Planck is also shown 
(black dashed). 



Figure 9: Varying f under n s = 1.04. a^/nf 
is fixed to 0.45 (blue), 0.40 (green), and 0.10 
(red), corresponding to a ~ 0.03, 0.008, and 
8 x 10~ 5 , respectively. 



4 Self-Interacting Curvatons 



As a simple example of curvaton potentials that steepen more rapidly than a quadratic, 
in this section we explore curvatons possessing a mass term and an additional higher-order 
polynomial term [18-24]: 



V{&) = A 4 



(4.1) 
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where A and / are positive constants with mass dimension, and m is an even integer with 
m > 2. / denotes where the higher-order term becomes important, while A sets the overall 
scale of the curvaton potential. Furthermore, the curvaton mass at the potential minimum 
is given by M 2 = 2A 4 // 2 . We especially focus on the region \a\ < /, where the analytic 
expressions give good estimations of the resulting curvature perturbations and a large running 
of /nl can show up. The validity of the analytic expressions at large a values are discussed 
towards the end of this section. 



The potential (4.1) gives 



4 A 4 J m(m — 1) fa* 



m-2 




^- l -3H!r-{ l + ^-^{j ] >• (42) 

_ = _ 4 m(m - l ){m -2) _V_(^ ] {l + .^ [ ._. ] h (4 . 3) 

where one can see that the resulting n s and a are positive and negative, respectively, for 
er* 7^ 0. We also note that — ,~ a ,, 2 takes its maximum value ^^ at (^f) m ~ 2 = , 2 .. . 

* ' {n s — iy 8 \ / ' m(m— 3) 

Hence especially for m < 10, then —a < (n s — l) 2 . 

Let us now examine cr osc since understanding its behavior in terms of a* is essential for 
discussing largely scale-dependent /nl produced from self-interacting curvatons. The Hubble 
parameter at the onset of oscillation (2.2) is 

1 A 4 / „ / a a 



H ^ = 1 c p[ 2 + m [T ] "' M!i 



thus the relation (2.22) between a* and a osc becomes 



In 



m + 2(^//)- (m - 2) 
m + 2(a osc /fy {m - 2) 



2., A 4 1 / /<7„ ^ m ' 2 



-(•n-2){- 3 K J ^ f2+ — 3 [2 + myf 



-1 



(4.5) 

By solving this equation, a osc is obtained as a function of a*. (Here a osc should take the same 
sign as a*, cf. Footnote #6.) We do this numerically, as shown in Figure 11, but the behavior 
of a osc can simply be understood by rewriting (4.5) as* 13 

2\°-j\ ^ ^ = -m+(m + 2(^) '" ^ | x 



cxp 



A 4 1 / (a. 



■("•-2) oM^ + n 2 



3 *tf 2 lf / 2 c-3 V V / 



osc 

m " 



m-2 x 1 



• (4-6) 



# 13 Here we exponentiate both sides of the equation in order to make clear the behavior of its solution, but 
let us note that since (2.22) itself is an approximate relation, one should in general be careful about the size 
of the errors when exponentiating the equation. 
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Since the left hand side of this equation is positive, a osc should not take values that make 
the right hand side negative, which can happen for e.g., large m, |cr osc //|, and A/"*A 4 /if^ f / 2 . 
Such values of <7 OS c are displayed in Figure 10, where regions of the cr osc / f - m plane on the 
right sides of the lines give negative values to the right hand side of (4.6). For example, 
when m = 8 and A/"*A 4 / 'Hf Q J 2 = 1, then a osc > 0.3/ cannot be a solution of (4.6), i.e. the 
curvaton rolls down before it starts its oscillations to values smaller than about 0.3/ even if a 
during inflation takes much larger field values* 14 . Such behavior is in contrast to a quadratic 
curvaton, for which one can check that cr osc oc cr*. In Figure 11 we plot a osc as a function of 
<r* by numerically solving (4.5) for m — 8 and A/"*A 4 /iJ^ f / 2 = 1. One sees that the growing 
rate of <t osc is suppressed at a* ~ / where the curvaton potential steepens due to the higher- 
order self-interaction, and that cr osc is limited to values smaller than ~ 0.3/ as is indicated 
in Figure 10. Such flattening of a osc is more significant for larger values of m and A 4 /H? nf f 2 
which make the potential steeper, and for larger N* providing a longer period for the curvaton 
to roll down during inflation. 
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Figure 11: a osc as a function of a*, for m = 8 
Figure 10: The right hand side of (4.6) and ^^4/^ /2 = L Solid and dashed lines 

becomes negative on the right sides of , , n /r> , - , . , 

° ° denote c = 9/2 and 5, respectively 

the lines. Blue and red lines denote 
ACAV-Htaf/ 2 = 1 and 10- 2 , while solid 
and dashed lines denote c = 9/2 and 5, 
respectively. 

The rolling of the curvaton due to the steep potential makes <t osc substantially smaller 



#14 Of course, discussions here are basically limited to curvaton field values that satisfy the approximation 
(2.1) until the onset of the oscillations. 
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than er*, and as we will soon see, this leads to a running non-Gaussianity. Let us now show 
the density perturbations generated by a self-interacting curvaton with the parameter set 
m = 8, A = 2.5 x 10 12 GeV, / = 4.3 x 10 13 GeV, and T a = 1(T 15 GeV (i.e. the energy 
density of the universe at the curvaton decay is pj ec ~ 65 GeV), along with the inflationary 
parameters H m { = 10 12 GeV, T^ = 3.8 x 10 3 GeV (i.e. the energy density at reheating (= 
inflaton decay) is p vch ~ 1.3 x 10 11 GeV), and jV* = 50 (see also Footnote #9). This set of 
parameters gives A 4 /if 2 lf / 2 w 0.02 (which determines the magnitude of n s — 1), and realizes 
the COBE (WMAP) normalization value as well as f ~ 10 -2 at around a*/f ~ 0.6, where 
rif NL blows up. The curvaton starts its oscillation before reheating, thus c = 9/2. The value 
of <7 osc in terms of <r* for this parameter set is illustrated by the blue solid line in Figure 11. 
Along with the analytic calculations, we also carried out numerical computations for the 
density perturbations in a similar fashion as we have done for NG curvatons in Section 3. 
The results are shown as functions of cr*/f in Figures 12 - 18, where the blue solid lines 
denote the analytic calculations and the blue dots the numerically computed results. Upon 
the numerical computations we have fixed the inflaton energy density to a constant during 
inflation, thus plotted n s and a, and the numerically computed n s and a in the same figures. 
In the presence of nonzero H* or H*, the actual spectral index and its running obtain offsets 
as shown in (2.26) and (2.27). For example, a 6 type chaotic inflation with J\f* = 50 realizes 
n s = n s — 0.06 #15 , and thus red-tilts the resulting perturbation spectrum at cr*/f ~ 0.6 where 
rif NL blows up. 

We also note that for the above parameter set, the curvaton 's classical rolling in a Hubble 
time is larger than its quantum fluctuations during inflation, i.e. V'/3H? ni > H in f/27r, at 
around a*/ f ~ 0.6 and larger. Although this is not necessarily the case for smaller a* (for 
a */f ~ 0-54 the quantum fluctuations become dominant over the classical rolling by the 
end of inflation), we have plotted in the figures down to small o* regions supposing that 
the curvaton dynamics during inflation can be treated as classical, in order to see the full 
(classical) ^-dependence of the density perturbations* 16 . 

The cr* <C / region is well-approximated by the familiar quadratic curvaton, but as one 
goes towards larger <r*, the system starts to behave quite differently. The power spectrum of 
the linear order perturbations (2.5) in Figure 15 increases in the small field regime with a* 
due to the increase of f. However for larger field values, the power spectrum starts to de- 
crease as a function of cr*. This is attributed to the behavior of a osc shown in Figure 11: 
When a*// approaches unity, cr osc ceases to grow as rapidly as a*, hence the terms in the 
expression (2.6) are insensitive to the value of a*, except through the explicit dependence 



# 15 Although the </> 6 chaotic inflation model is excluded due to too large tensor-to-scalar ratio and red-tilted 
spectral index when the inflaton is totally responsible for the primordial density perturbation, in the curvaton 
framework, chaotic inflation models with high order polynomials are still viable. 

# 16 We also remark that for the example parameters we have adopted, the classical rolling is not so large 
compared to the quantum fluctuations at around er*// sa 0.6, hence random effects from quantum fluctuations 
may also need to be taken into account for a more rigorous treatment. 



dN /da* oc da osc /da* oc 1/V'(cr*). As a consequence, the power spectrum becomes a decreas- 
ing function of a* for large field regime. In other words, in the large a* regime the steep 
potential forces the curvaton to roll down to take similar values for cr osc almost independently 
of a*, therefore dj\f*/da* is suppressed. The peak of the power spectrum is located in the 
intermediate region, i.e. at a* ~ 0.6/, and this roughly matches with the asymptotic value 
of cr osc that can be read off from Figure 11 (or 10). When starting from large a*, the initial 
field fluctuations are reduced as the curvaton rolls down to field values below the asymptotic 
value of dose, and thus the linear perturbations are suppressed. 

The existence of the peak in the power spectrum indicates that the non-linearity parameter 
/nl crosses zero and its running n/ NL blows up. This is understood also in terms of f\ and 
/ 2 (2.15), which are plotted in addition to /nl m Figure 16. At o* <C f the amplitudes of 
/i and / 2 are comparable and their sum results in a positive /nl as is the case for quadratic 
curvatons (cf. Footnote #5). However as one takes larger a*, / 2 which is a function of <t osc and 
f becomes insensitive to a*, while /i rapidly decreases on the negative side through its explicit 
dependence on a*, resulting in a negative /nl- One clearly sees that the condition (2.35) is 
satisfied at around a*// ~ 0.6 where /nl crosses zero, being able to produce a large n/ NL even 
under a suppressed a. Furthermore, since f<l, the amplitude of /nl itself is large (except 
for the very vicinity of its vanishing point) while satisfying (2.35), i.e. 1 <C |/nl| "C |/i|. 
(This is in contrast to quadratic or NG curvatons, for which /nl vanishes at rather large f. 
This was why the product nj NL / NL from NG curvatons was suppressed even when n/ NL blew 
up, cf. Figure 9.) Therefore a large and strongly scale-dependent /nl can be produced, whose 
running is in the detectable range by upcoming CMB observations. Contour lines in the rif NL 
- /nl plane are shown in Figure 18, where the black dashed line denotes the Planck detection 
limit ^/ nl /nl/50 = 0.10, cf. (2.36). Compared to the Figures 12 - 17, in Figure 18 we have 
added more numerically computed points, which are equally spaced in terms of Act*. 

Let us also comment on the validity of the analytic estimations. In most of the displayed 
region, the results from the analytic calculations are in good agreement with those from the 
numerical computations. However, as one goes towards larger a*, i.e. larger a osc , a period 
of non-sinusoidal oscillations along the a m potential needs to be taken into account. Thus 
the deviations between the analytic and numerical results which one can already see in, 
e.g. Figure 16, becomes even larger for ct* beyond the region displayed in the figures. In 
Appendix A we extend the expressions in Section 2 to incorporate a period of non-sinusoidal 
oscillations, which becomes important for large cr osc values. We also note that, with the choice 
of the parameter set here, the curvaton's effective mass during inflation is comparable to H in f 
for a* ~ /, which violates the slow-roll approximation and sources additional errors to the 
analytic estimations. 

For self-interacting curvatons with even larger cr osc , we should remark that our approx- 
imations for the analytic expressions can break down. This is because, especially for cases 
with large m and cr osc //, the potential curvature V" quickly decreases after oscillation starts 
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(which we have defined by when \&/Ha\ = 1) as the curvaton rolls down to smaller field val- 
ues. This forces the curvaton to recover (though not completely) the attractor dynamics (2.1) 
for a short time after t osc . As a consequence, the simple picture of the curvaton suddenly 
starting its oscillations at t — t osc is no longer a good approximation. Such breakdown of the 
approximations happen especially for large m or cr*/f, or small A 4 /H 2 nf f 2 , which realize large 
values for <r osc . We also note that large m and a* can lead to the breakdown of the attrac- 
tor solution (2.1) well before the onset of the oscillations, which can also spoil the analytic 
estimations. Density perturbations from a self-interacting curvaton with large a* values have 
been worked out numerically in [33]. 

In Figures 19 and 20 we further show contours in the n^ NL - /nl planes for a varying a* 
under fixed values of m, A 4 /H 2 nf f 2 , and f (instead of fixing individual parameters). We have 
taken A/"* = 50 and t osc < ^reh, i-e. c = 9/2, but we remark that the results do not change 
significantly for t osc > £ re h, i.e. c = 5. (It should also be noted that in order to compute /nl 
and nf NL , we only need to specify m, A 4 /H 2 n{ f 2 , f, a*/ f, A/"*, and c. Note also that in order 
to determine the liner perturbation amplitude V^, we further need to fix the ratio A//.) Here 
we are fixing f, which can be considered as varying r CT along with cr*/f, hence the resulting 
/nl and rif NL can behave somewhat differently compared to Figures 16 - 18, especially in the 
region |cr*| <C /. 

Steeper potentials, i.e., larger A 4 /H 2 ni f 2 (or larger m), realize flatter functions of <r osc in 
terms of a*, which tend to produce /nl with stronger scale-dependence as shown in Figure 19. 
Moreover, smaller f basically leads to larger |/nl|, pushing rif Nh towards the detectable regions 
as one sees in Figure 20. A somewhat different behavior is seen for the f 3> 1 case shown as the 
red contour line in Figure 20, where /nl is always negative and the <r* — > limit corresponds 
to the endpoint of the line at nj Nh = 0. This is understood as the f ^> 1 behavior of quadratic 
curvatons, cf. (2.19). Moreover, we have plotted the red line up to cr*/f ~ 0.9, beyond 
which n s becomes larger than 1.5 and the curvaton no longer slow-rolls during inflation. This 
corresponds to the other endpoint of the red line, where rather large \nf NL \ (together with 
a largely negative /nl) is realized due the large h s — 1 and \a\, coming close to satisfy the 
condition (2.34). 

In summary, the steep potential forces the self-interacting curvaton to roll down to small 
field values by the onset of the curvaton oscillation. This greatly diminishes the initial dif- 
ferences in the curvaton field values during inflation a*, thus suppresses the resulting linear 
order density perturbations. When fixing all the parameters of the system except for a*, 
then a maximally large linear perturbation amplitude is obtained from a* that is close to the 
asymptotic value of a osc , and around this a* value is where a strongly scale-dependent /nl 
is produced. We expect such behavior of a self-interacting curvaton to be a rather generic 
feature of curvatons whose potentials are approximated by quadratic around their minimum, 
but steepens more rapidly than a quadratic away from the minimum. 

20 



a 



1.15 



1.10 



1.05 



1.00 




0.2 0.4 0.6 0.8 / 



0.000 
-0.005 
-0.010 
-0.015 
-0.020 ; 
-0.025 : 
-0.030 



0.2 




a-,,, 



f 



Figure 12: Curvaton contribution to the 
spectral index (2.28). 




Figure 13: Curvaton contribution to the run- 
ning of the spectral index (2.29). 




Figure 14: Energy fraction at decay. 

1000 r 



500 



-500 

1000 



Figure 15: Linear order perturbations. 
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Figure 16: Non-Gaussianity. /nl ; blue solid, 
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Figure 17: Running of non-Gaussianity. 
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Figure 18: Plot in the n/ NL - /nl plane. The region 0.45 < cr*/f < 0.65 is shown in the 
displayed area. The expected observational sensitivity of Planck is also shown (black dashed). 
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Figure 19: Varying ajf for A 4 /# 2 nf / 2 = Figure 20: Varying ajf for f = 10 3 (red), 
0.03 (red), 0.02 (blue), 0.01 (green), with 1(T 2 (blue), 10" 22 (green), with fixed m = 8 
fixed m = 8 and f = 10~ 2 . The expected ob- and A 4 /if 2 nf / 2 = 0.02. 
servational sensitivity of Planck is also shown 
(black dashed). 
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5 Mixed Curvaton and Inflaton 

In the preceding two sections, we have assumed that fluctuations from the inflaton are neg- 
ligible and those from the curvaton are only responsible for cosmic density fluctuations. But 
now in this section, we further consider the case where the inflaton also contributes to the 
density perturbations as well as those from the curvaton, which is called "mixed curvaton and 
inflaton model" in the literature [34-38]. 

We start with giving the expression of the curvature perturbations in terms of the SJ\f- 
formalism, assuming that the field fluctuations of the inflaton 50 and the curvaton 5a obey 
Gaussian distributions (to be more precise, in the sense discussed below (2.13)), with Vs^ik) = 
Vsa{k) = (if |fc= a _ff/27r) 2 when k exits the horizon. We further suppose that there are no direct 
couplings between the curvaton and inflation, thus no correlations between their fluctuations 
5(j) and 5a. Then, considering a homogeneous and isotropic universe background, the power 
spectrum Vq and the non-linearity parameter /nl at the CMB scale in this model can be 
written as [37] 

P c =(A/J+A/; 2 )(§) 2 , (5.1) 

5 A/; 2 A/^ + Af^Ka + 2M a M^ a 

at the leading order in terms of the field fluctuations. Here the subscripts and a denote 
partial derivatives with respect to the fields, i.e. Af<j, = dj\f/d(j)*, Af a = dAf/da*, etc. 

Given that the end of inflation is independent of 0* (so that quantities such as the energy 
density of the universe and the inflaton field value at the end of inflation are independent 
of 0*), one can find that 

M+ = -X 1 . ( 5 - 3 ) 

0* 

where we have assumed that H and during inflation are determined merely by (e.g., as 
in slow-roll inflation). Here we should remark that there are actually further contributions to 
Mcj> since fluctuations in the duration of inflation due to 5(f) lead to curvaton field fluctuations 
in the post-inflationary era, which can be seen as 

da osc i V {a osc ) 

— — = (1 - X{a osc )) . . 5.4 

•90* 3#*0* 

This can be derived in a similar fashion as (2.23), considering a curvaton potential V(a) that 
is a function only of a. In this sense, the inflaton can further produce density perturbations 
after inflation by fluctuating the curvaton. Moreover, we note that such effect sources Af^, 
which otherwise vanishes since the inflaton and curvaton only affect the e-folding numbers 
during and after inflation, respectively. However, in this section, we ignore effects due to (5.4) 
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since we would like to show that /nl can become largely scale-dependent simply by having 
multiple sources (even in the absence of the cross term A/^), and also because their effects 
are negligible for the case of a quadratic curvaton and a slow-rolling inflaton which we will 
study later in this section* 17 . 

For convenience, we introduce a parameter which represents the fraction of the curvaton 
contribution to the total power spectrum: 



m 



m+N? 



(5.5) 



This quantity satisfies < q < 1 by definition. Then using this quantity, the power spectrum 
and the non-linearity parameter /nl can be rewritten as 
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Now let us write down the expressions in terms of the effective potentials. We assume that 
the inflaton and the curvaton are decoupled from each other so that the total potential can 
be divided into the inflaton potential U(<p) which drives inflation as well as determines the 
inflaton dynamics, and the curvaton potential V(a) which governs the curvaton dynamics. 
Both potentials are considered to have no explicit dependence on time. Furthermore, for the 
inflaton, we assume a slow-roll inflation, i.e. 

3H<j) = -U', 3H 2 M* = U, (5.8) 

where a prime on U denotes differentiation with respect to <fi. These approximations give 
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Introducing the slow-roll parameters 
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(5.9) 
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^ 17 To be precise, for the mixed case of a quadratic curvaton (with mass m) and a slow-rolling inflaton, effects 
from (5.4) on the power spectrum Vq are suppressed by a factor m 2 /H% compared to the contribution (5.3). 
As for /nl, unless f takes a specific value f rj 1.85 which vanishes MaajM 2 , contributions from (5.4) on /nl 
are either sufficiently smaller than unity, or subleading compared to other terms. 
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The curvaton is also assumed to slow-roll during inflation, i.e. 3H& = —V, where a prime 
here for V denotes a derivative in terms of a. For the 5M from the curvaton we adopt the 
formulae used in Section 2: N a is shown in (2.6), and M aa can be read off from (2.14). 

Then, at the leading order in terms of the slow-roll parameters e* and 77*, we obtain 
expressions for the spectral index n s — 1, its running a, and the product ?i./ nl ,/nl as follows: 



n s - 1 ~ (1 - g)(-6e* + 2^) + q -2e* + 2 



a ~(1 - g) (-24e 2 + 16e*r?* - 2£*) 



, V"M 



(5.12) 



+ 4g(l - q) 2e* - 77* + 
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(5.14) 
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g 2 r"(<7* 



AC 3^2 
where we have further introduced another slow-roll parameter for the inflaton: 

U'U'" 



* - M P -[P 



(5.15) 



Let us stress that upon obtaining the above results, we have made use of the slow-roll 
approximations for and a, and further used the analytic formulae for M a and N aa which 
also relies on some approximations (including the sudden decay) as discussed in Section 2. 
It should be noted that errors contained in the approximations can accumulate, especially as 
one goes to higher order derivatives. 

The running of /nl receives additional contributions in the mixed case from the second 
line of (5.14), which is absent for a pure curvaton or inflaton case. In particular, large /nl and 
rif NL can be obtained when the curvaton mainly contributes to the (running) non-Gaussianity, 
while the inflaton mainly generates the liner order perturbations. In fact, the running a also 
has the contribution specific to the mixed model such as the second line of (5.13). However, it 
is suppressed as (0(e*,?7*,77 ff )) 2 with r\ a = V /(3H%) corresponding to a slow-roll parameter 
for a, hence a is small in this model. 

To see to what extent the running of non-Gaussianity can be large in this kind of model, 
as an example, let us consider a quadratic curvaton V{a) = ^m 2 a 2 for which nj NL vanishes 
for the curvaton-only case, i.e. q = 1. (Note that for quadratic curvatons, /nl and n/ NL are 
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Figure 21: Contours for f = 10~ 2 (red line) 
and 10 -3 (green line) with q being varied in 
the n/ NL -/NL plane. Here we assume that 
m 2 /(3H 2 ) = 0.05. The slow-roll parame- 
ters for the inflaton sector are taken as e* = 
0.016, ?7* = 0.025, and £„ = 0.0004. The ex- 
pected observational sensitivity of Planck is 
also shown (black dashed). 



Figure 22: The same cases as Figure 21 are 
shown as a function of q. /nl is plotted for 
f = 10~ 2 (red) and 10 -3 (green). ri/ NL is 
for f = 10~ 2 (red dashed) and f = 10~ 3 
(green dashed). Notice that the value of n/ NL 
becomes insensitive to f as q approaches to 
unity. 



independent of a* and <J osc .) By doing this, we can see n^ NL generated from the multi-source 
nature of this scenario. 

Given that e*, |r/*|, |£*|, m 2 /H 2 <C 1, then one can check that |/nl| ~^> 1 along with 
I^/nl/nlI ^ 1 are realized if and only if 



<9 2 (l-<?) 



2e* - 77* + 
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3H* 



(5.16) 



is satisfied. In such a case, /nl and its running are approximated by 
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(5.17) 



In Figure 21, plots in the 71/ nl -/nl plane are shown for cases with f = 10 -2 and 10~ 3 
with q being varied. For the curvaton mass, we assume m 2 /(3if 2 ) = 0.05. The slow-roll 
parameters for the inflaton sector are taken as e* = 0.016, 77* = 0.025, and £* = 0.0004, which 
correspond to the values for 4 chaotic inflation model with TV* = 60. Since non-Gaussianity 
mainly comes from the curvaton sector, /nl becomes larger as q approaches unity as read off 
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Figure 23: Contours of /nl and % NL in 
a^/Mp-r plane. Here we show the contours 
of rif NL = 0.2 (green) and 0.1 (blue). The 
values of /nl are shown in the figure. Other 
parameters are assumed as m 2 /(3H%) = 0.05, 
e* = 0.016,7/, = 0.025, and & = 0.0004. 



Figure 24: Shown are contours of n,f NL f NL = 
5 (green) and 2.5 (blue), which correspond 
to the detection limit values for Planck and 
CMBPol, respectively. Other parameters are 
taken as the same as in Figure 23. 



from (5.17). But on the other hand, nf NL becomes small in this limit because the running 
of the non-Gaussianity arises due to the multi-source nature of the model and the case with 
q — y 1 corresponds to the single- source limit. 

Regarding the scale dependence of /nl, n f Nh takes its maximum value at intermediate 
range of q since one has to have non-zero value of g(l — q) to have sizable n^ NL in this kind of 
a mixed model. It should be noted here that although the value of /nl becomes very small 
as q goes to 0, however, if we assume a very small value f by choosing some appropriate 
parameters for the curvaton sector, the smallness of q is somewhat compensated and /nl can 
be sizable. Thus for some values of q, both / NL and nj NL can be large enough to be detected 
by Planck satellite as seen from Figure 21. To see how /nl and rif NL depend on q, we also 
plot the values of these quantities as a function of q in Figure 22, from which we can clearly 
see how the size of q affects these quantities. 

Since q is determined once we fix e#, a, and f, we also show the prediction of the model in 
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the o*-r plane in Figure 23 where contours of /nl and rif NL are shown. As already mentioned 
in the previous section, the detectability of rif NL rather depends on the product of n,f NL f NL , 
hence we further plot contours of this quantity in Figure 24. In the figure, the contours 
corresponding to the detection sensitivity limit for Planck and CMBPol are depicted. From 
the figure, we can see that, in some parameter region, rif NL can be detectable in near future 
observations. 

6 Multi-curvaton 

As a final example model related to the curvaton, we in this section consider the multi- 
curvaton scenario where multiple curvaton fields can be responsible for cosmic density per- 
turbations. Models of this kind have been investigated in the literature [39, 40] and some 
explicit examples in particle physics have also been discussed [52]. In this section, we again 
neglect the contribution from the inflaton fluctuations to the total curvature perturbations. 
Although there may generally exist many curvaton fields, here we consider a two curvaton 
case. 

As in the previous sections, first we write down the expressions for the spectral index n s , 
the running a, non-linearity parameter /nl 5 and its running rif NL . In general, the curvature 
perturbations in the model can be written as, up to the second order of the field fluctuations, 

C = M a 8a* + M b 8K + - {M aa (8a*) 2 + M bb (8b*) 2 + 2Af ab 8aJh) , (6.1) 

where 8a* and 8b* are fluctuations of the curvaton fields "a" and "6" , respectively, and M a = 
dM /da* and so on. Throughout this section, we denote "a" and "6" curvatons as the ones 
that decay first and later, respectively. 

The explicit forms of M a , M b , etc., can be expressed in terms of the curvaton potentials 
by extending the discussions in Section 2 to the multiple curvaton case. However, since 
the explicit expressions are in general very complicated and we are interested in the scale- 
dependence which comes from the multi-source nature of this scenario, here we limit ourselves 
to the case where the curvaton potentials are quadratic, i.e., 

V(a,b)= 1 -my+ 1 -mlb 2 . (6.2) 

We do not specify which of the masses m a and m b is larger. We further suppose that the 
curvatons start to oscillate after the inflaton decay, and that the energy density of the cur- 
vatons are negligibly tiny compared to the total energy of the universe until both curvatons 
start their oscillations. Then, one can find that M a , M aa , • • • are of the form 

M a oc — , M aa oc — , M ab oc — — , (6.3) 

(the same for replacing a with b), where the constants of proportionality are given as functions 
of the ratios between the energy densities of a, b, and radiation upon decays of a and b (the 
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full expressions of the coefficients up to the third order are given in [38], where the density 
perturbations are computed in a different approach from Section 2): 
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(6.4) 
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Here we have defined f a i and f b2 similarly to f defined in (2.8). p a , Pb and p r are energy 
densities of a curvaton, b curvaton and radiation, respectively, while t a d cc and t b d ec are the 
times at the decay of a and b curvatons. Notice that f al and f b2 are defined at different times, 
thus we also put the subscripts "1" and "2" in addition to "a" and "6" in the definitions. 
We also remark that the p r in the denominators of f al and r b i come from the inflaton decay, 
while p r in fb2 also includes the decay products from the a curvaton. 

Since the proportionality factors in (6.3) are determined by the three ratios (6.4), the scale- 
dependence of the density perturbations are sourced through a* and &*. Hence the spectral 
index n s and its running a for the power spectrum are (assuming Gaussian field fluctuations 
in the sense explained below (2.13), with Vs a (k) = Vsb(k) = (H\k= a H /2ti) 2 when k exits the 
horizon, and no correlations between 8a and 8b) 
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TT 

n s - 1 ~ 2 -J + 2K a r] a + 2K b r] b , 



2 fl " 4 fi " A Jf2 ( K ^a + K bVb ) + ±K a K b ( Va - Vb ) 2 



(6.5) 
(6.6) 



where r\ a and r\ b are defined as 
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We have also introduced the parameters K a and K b representing the fractional contribution 
of the curvatons a and b to the total power spectrum: 
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which clearly satisfies K a + K b = 1. 

Regarding the non-Gaussianity, the non-linearity parameter /nl is generally given by 
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and its running n/ NL can be computed as 
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(6.10) 
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In the limit of either one of these curvatons being solely responsible for density fluctuations 
(i.e., in the limit of K a (K b ) — > 1 and K b (K a ) — > 0), the above formulae reduce to a single 
curvaton case. What is specific to this model is the terms with K a K b which disappear in a 
single curvaton scenario. In particular, when the masses of two curvatons are different (i.e., 
Va 7^ Vb), the running rif NL can be sizable, which has been pointed out in [29,30]. In fact, 
a also receives the contribution which originates from the mass difference of the curvatons. 
However, as seen from (6.6), such a term is suppressed with (rj a — r] b ) 2 , thus a would be small 
in this model. 

In fact, even if we assume a quadratic potential for the curvatons, general expressions for 
the coefficients such as M a) M aa and so on are still very complicated. Thus in the following 
we discuss some limiting cases where both curvatons are subdominant and dominant at the 
time of their decays, paying particular attention to the non-linearity parameter fah and its 
running n/ NL . 

6.1 Both curvaton subdominant at their decays 

First we consider the case where the energy densities of both curvatons at their decays are 
subdominant, i.e., 

f a i, r m < 1. (6.11) 

Since we adopt a purely quadratic potential for the curvatons, we can treat the curvatons as 
nonrelativistic fluids once they start to oscillate. By using the sudden decay approximation, 
one can show that the explicit forms of (6.3) in this case become [38] * 18 

M fal l AT ^ l AT fal X AT ^ al f b2 1 f b2 1 

iV a ~- , N b ~ — — , N aa ~—- — , N ab ~ -, N bb ~— — 6.12 

2 a* 2 6* 2 a% 4 ajo* 2 b% 

which give n s ,a,/NL, and nf NL through the formulae given in (6.5)-(6.10). For convenience, 
we also introduce a parameter representing the ratio of the amplitude of the a and b curvatons: 

A = — . (6.13) 

a* 

In Figure 25, we plot /nl and n,f NL as functions of A for several parameter sets of (f al , 
^6i, ^62) Va, Vb)* 19 - Actually, these parameters are not independent, as will be discussed 
around (6.16). Here we take different values for r\ a and rj b ^ 20 , which can give a sizable value 



# 18 To be precise, here we are also assuming that the decays of the curvatons a and b are well separated along 
the expansion timeline. 

#19 p or completeness, we also give the value of fbi used in plotting the figures since we use a full expression 
provided in [38] for the computations. However, we remark that its actual value is not important in this limit 
as can be noticed from (6.12) where an explicit dependence on f b i disappears. 

# 20 As can be seen from (6.16), the parameters are not independent and r\ b in the figure takes values much 
smaller than r\ a . However, when i] b -C rj a is satisfied, the explicit value of rj b is not important for obtaining 
/nl and n /NL . 
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for rif NL as noted above. When fixing f al and f^, the fraction K a increases together with 
A. From (6.10), one can see that large values of n/ NL /NL would be obtained when K a ~ 
Kb (corresponding to A ~ r^/fai ~ O(10 2 ) for the parameter sets assumed in Figure 25), 
otherwise the factor K a Kb becomes very small. Furthermore, in the limit of both curvatons 
subdominant we are considering here, one has 

Ka J?_ Mbb_ _2_ Kb , , 

AC 2 "fal' A/?~f 62 ' KM b ~ [ } 

Thus, when f a i € f(, 2 « 1 and K a ~ i^ where both /nl and n/ NL can be large, the 
terms originating solely from a curvaton dominate in (6.9) and (6.10) (i.e., K^(M aa /Af 2 ) and 
K a {N a a/Na) in the expressions of /nl and nj Nh , respectively). In this case, we can make a 
rough estimate of the values of /nl and n/ NL as 

/nl ^ |^ a 2 ^, n /NL ~ 4(1 - K a )( Va - rj b ). (6.15) 

<J ' al 

In fact, for the case with (f al ,f5 2 ) — (10~ 4 , 10~ 2 ), rif NL takes its maximum value at around 
K a ~ 0.1 (fCf, ~ 0.9). Since now we are assuming that f a \ <^ 1, large /nl is possible even 
with a relatively small value of K a . Furthermore, although rj a gives a positive contribution to 
n s (see (6.5)), small i^ a together with rji ~ imply that the spectral index can be red-tilted 
when assuming a large field inflation model with H*/H 2 ~ —0(0.01). Thus, in this model, 
we can have large /nl and rif NL with red-tilted power spectrum for some parameter region. 

As emphasized several times in the previous sections, the detectability of n/ NL also depends 
on the size of /nl, thus we also show a plot in the n/ NL -/NL plane in Figure 26. For reference, 
the detection limit lines for Planck |^/ nl /nl| = 5 are also shown. In the figure, the value of 
A is varied as in Figure 25. In some region, the predictions of the model for / NL and rif NL 
are beyond the detectable sensitivity line, thus if ri/ NL is detected in future experiments, this 
model may also be a target of serious study. 

We should also remark that in the case discussed above, a detectable level of rif NL can 
be obtained at the expense of fine-tuning of the parameters. Since here we are assuming 
that both curvatons begin their oscillations during the radiation dominated epoch after the 
inflaton reheating, the ratio of the energy densities between a and b curvatons at the time of 
a curvaton decay is given by (independently of whether m a ^ m&) 



?"al _ Pa 

hi Pb 






where for simplicity we have taken b osc /a osc ~ 6* /a*. Therefore one sees that the parameter 
sets chosen in Figure 25 suppose the mass m b to be extremely suppressed, much smaller 
than m a . See also discussions in Footnote #20. 
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Figure 25: Plots of /nl (left panel) and rif NL (right panel) as a function of A for the cases 
with both curvatons being subdominant at their decays. Here we fix the parameters as 
(r a i,f b i,r b 2,Va) = (l(T 4 ,iO- 5 , 0.07, 0.05) (red) and (10" 3 , 10" 4 , 0.5, 0.02) (green), while rj b 
takes values much smaller than i] a (cf. Footnote #20). 
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Figure 26: The same cases as Figure 25 are plotted in the /nl _ ^/ nl plane. The expected 
observational sensitivity of Planck is also shown (black dashed). 
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6.2 Both curvaton dominant at their decays 

Now we consider the case where both curvatons are dominant at the time of their decays, i.e., 

f a i, r b 2 > 1, f i > hi- (6.17) 

In this case, one obtains 
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It should be noted here that £ depends on f& 2 but not on f al nor f&i at leading order as far 
as we assume the condition (6.17). Also one can easily see that 
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(6.19) 



Hence, when the assumption of (6.17) holds and K a ~ K b where nj NL can be sizable, /nl and 
?2/ NL are roughly estimated as 
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(6.20) 
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Figure 27: Plots of / NL (left panel) and rif NL (right panel) as a function of A for the 
cases with both curvatons being dominant at their decays. Here we fix the parameters as 
(r a i,f b hh2,Va) = (10 3 ,10 2 ,10 3 ,0.05) (red) and (10 3 ,10 2 ,2x 10 3 ,0.02) (green), while rj b takes 
values much smaller than r\ a (cf. Footnote #20). 
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Figure 28: The same cases as Figure 27 are plotted in the /nl _ ?V nl plane. The expected 
observational sensitivity of Planck is also shown (black dashed). 



In Figure 27, we plot /nl and rif NL as functions of A for several parameter sets of (f al , f bl , 
^62, Vai Vb)- As mentioned above, as long as we assume (6.17), the expressions for Af a ,Af aa 
and so on are given as (6.18). Thus the assumption on the values of r a \ and r&i have little 
effects on /nl and n/ NL . Under a fixed f& 2 , the fraction K a increases along with A. It should 
be noticed that /nl becomes negative when A is small, i.e. K a <C 1, due to the negative 
values of A/ftb/A^ 2 . However, as A becomes larger (K a is larger), / NL crosses zero and increases 
to large positive values as seen from Figure 27. Hence nj NL blows up at around the region 
where /nl ~ similarly to the case of the self-interacting curvaton. Further increasing A, 
then K a becomes comparable to Kb and (6.20) is satisfied. For even larger A, n/ NL again 
becomes small. To see the detectability of rif NL in this case, we again show the plot in the 
tij nl -/nl plane in Figure 28, from which we can clearly see that rif Nh can be detectable in 
some parameter region as in the subdominant case discussed in the previous subsection. We 
remark that n s can also be red-tilted at around the region where n/ NL becomes large. In 
fact, with the parameter sets assumed in Figures 27 and 28, the contribution to the power 
spectrum from the a curvaton is somewhat small as K a ~ 0(0.01) when n/ NL blows up. There 
the positive contribution of K a i] a in n s becomes small, in which the spectral index can be 
almost given by n s — 1 ~ IH^jHl and be red-tilted for large field inflation models. 

We also note that (6.16) requires the mass m\, to be highly suppressed in the above case, 
as in the previous subsection. Here, especially since the a curvaton is dominant at its decay, 
which means that p a > p^, the constraint (6.16) indicates \/m a /mb > A 2 . As seen from 
Figure 27, the product of /nl^/ nl can be large when A ~ (9(10 2 ), which requires a large 
hierarchy between m a and m?,. 
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7 Conclusions 

In this paper we have investigated the scale-dependence of linear and second order density 
perturbations produced in the curvaton mechanism. In particular, we focused on the runnings 
such as a and n/ NL . Since models of the curvaton have several variants, we discussed the 
representative models: curvatons with non-quadratic potentials, and multiple source cases. 
For the later model, we discussed multi-curvaton and mixed curvaton and inflaton ones. 

Non-quadratic curvaton potentials give scale-dependence to the non-Gaussianity, as well 
as to the spectral index of the linear perturbation spectrum. We have especially shown that 
the local-type /nl produced from curvatons can strongly depend on the scale, even when the 
linear order perturbations are nearly scale- invariant. We analytically computed the running 
n/ NL = d\a.\f^i\/d\n.k for curvatons with general energy potentials in (2.31), and obtained 
conditions (2.34) and (2.35) under which /nl can be strongly scale-dependent. 

The two conditions (2.34) and (2.35) can be satisfied respectively by curvaton potentials 
that flatten and steepen compared to a quadratic one. As an example of the former case, 
we looked into pseudo-Nambu-Goldstone curvatons with cosine-type potentials, for which it 
was shown that n/ NL is directly related to the amplitude of a and thus strictly constrained by 
current observational bounds on running spectral index, unless the inflationary mechanism 
realized a largely time-varying H such that cancelled out the contribution of a to the running 
spectral index. For the latter case, we studied self-interacting curvatons in detail and saw 
that the steep potential wipes out the initial field fluctuations 8a*, suppressing the resulting 
density perturbations. This lead to the production of a strongly scale-dependent /nl even 
for a suppressed running spectral index, when a* was away from, but not too away from the 
potential minimum. Moreover, the scale-dependence of /nl from self-interacting curvatons 
could be as large as to be detected by upcoming CMB experiments. 

In Sections 5 and 6, we further discussed the mixed curvaton and inflaton, and multi- 
curvaton models. In these models, due to the multiple source nature, large values of nj NL can 
also arise in a different way from the conditions discussed for the curvaton only case. We 
have made quantitative discussion for the running and shown in what cases n^ NL and /nl can 
be both large enough to be detected in future cosmological observations. 

The analytic methods used in this paper clarifies the intriguing behaviors of the density 
perturbations produced from non-quadratic and/or multiple curvatons, providing us with a 
systematic framework for studying the curvaton scenario in general. While for pseudo-Nambu- 
Goldstone curvatons a running / NL was accompanied by a large a and thus was constrained, 
the steep potential of self-interacting curvatons allowed strongly scale-dependent /nl- We 
expect that in terms of the /nl scale-dependence, these two examples respectively give typical 
behaviors of curvatons whose potentials flatten/steepen compared to a quadratic. Multiple 
source scenarios could also produce large n/ NL . The systematic approach we presented will 
be helpful for probing the physics of curvatons when combined with upcoming data. 



35 



Acknowledgements 

T.K. would like to thank Neal Dalai, Teruaki Suyama, Fuminobu Takahashi, and Jun'ichi 
Yokoyama for helpful conversations. The work of T.T. is supported in part by JSPS Grant- 
in- Aid for Scientific Research No. 23740195 and also by Saga University Dean's Grant 2011 
For Promising Young Researchers. 

A Density Perturbations from Curvatons with Non- 
Sinusoidal Oscillations 

The analytic expressions in Section 2 can be generalized to include a period of non-sinusoidal 
oscillations, as was discussed in Appendix B of [28]. Here we suppose that the energy density 
of the oscillating curvaton initially redshifts as p a oc a~ n with a constant n, then when 
its energy approaches a certain value p aS i n , the curvaton suddenly switches to a sinusoidal 
oscillation and then redshifts as p a oc a" 3 . We take p aS i n as a constant (i.e. independent 
of <7*), considering for e.g., self-interacting curvatons whose oscillations are determined by the 
curvaton field value. The curvaton energy density is assumed to be negligibly small until the 
inflaton decay or the onset of the curvaton oscillation, whichever is later. 
Then the linear perturbation amplitude becomes 



with 

M ^_ t f 3 Vja^) 3X(q osc ) | V'(a OBC ) 

da*-4 + 3r [i M ° osc)) \n V(a osc ) ' a osc J V'(<r.) ' l j 

where A(<t osc ) is defined in (2.9). The spectral index n s and its running a are the same as in 
(2.26) and (2.27), respectively. Furthermore, the non-linearity parameter is 
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(A.3) 



and its running takes the form 



-i 



1 5(4 + 3f) (3V'(a osc ) 3X(a osc ) \ l . w .. V'(a*) V 



r '"(o* 



n hL * 7^ v 1B " ' -T/T^T - -"^^ (1-*(*~) W V m ■ ( A - 4 ) 
/nl 18r {nV(a osc ) a osc J V'(a osc ) H* 
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Of course, the above equations reproduce (2.6), (2.14), and (2.31) for purely sinusoidal oscil- 
lations, i.e. n = 3. 

For example, for a self- interacting curvaton (4.1) with m — 4, then when a osc is larger 
than /, we need to use the expressions in this appendix with n = 4. 
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